Discrete symmetries in the Weyl expansion 
for quantum billiards 

0\ ■ Nicolas Pavloff 

Q ' Division de Physique TheoriqueQ, Institut de Physique Nucleaire, F-91406 Orsay Cedex, France 

in ■ 

abstract 
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f — ■ We consider 2 and 3 dimensional quantum billiards with discrete symmetries. The boundary 

^—1 . condition is either Dirichlet or Neumann. We derive the first terms of the Weyl expansion for the 
level density projected onto the irreducible representations of the symmetry group. The formulae 
require only the knowledge of the character table of the group and the geometrical properties (such 

■^ . as surface, perimeter etc ...) of sub-parts of the billiard invariant under a group transformation. As 

- ' an illustration, the method is applied to the icosahedral billiard. 
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Billiards have been extensively studied in the context of quantum chaos, both from the semi- 
classical point of view (periodic orbit theory, see e.g. Ref. |j|]) and using random matrix theory (for a 
review see pl). In these two approaches it is important to subtract the smooth part of the level density 
in order to study the oscillations (trace formula) or the short range fluctuations (random matrices) 
around the average behaviour. The asymptotic behaviour of the smooth part of the level density — 
characterized by the Weyl expansion — is also of interest in a great variety of other physical and 
mathematical problems (see Refs. @, §, |^, ^, 0, P)- 

In this paper we study two and three dimensional billiards : a quantum particle is enclosed in 
a compact region B of space and we impose to its wave-function %[) Dirichlet or Neumann boundary 
conditions on the frontier dB {dB may or may not be a smooth manifold) . Thus ^ verifies the following 
equation : 

(A + k'^)'ip = inside B , (1) 

^jJ = on dB (Dirichlet) 
n.'V'ip = on dB (Neumann) . 

We will assume that B remains invariant under a discrete point group Q. Then the eigenfunctions 
{ipn}n>o and eigenvalues {kn}n>o can be classified according to the different irreducible representations 
(irreps) of the group. We label these irreps by an index (a), each has a dimension d". Hence the total 
level density p{k) can be written as : 

pik) = J2p''{k), (2) 

a=l 

where r is the total number of non-equivalent irreps of Q, and p'^{k) is the density of levels belonging 
to irrep a : 

oo 

p-{k) = d" ^ 6ik - kZ) . (3) 

n=0 

p{k) (resp. p^ik)) can be separated in a smooth function of k, p{k) (resp. p"{k)) plus an oscillating 
part p{k) (resp. p"-{k)). The study of the oscillating part p°' needs a careful treatment and has been 
addressed in Ref. |9|, |l^, |l^, 12 1. Our aim is to write the first terms of an explicit Weyl expansion for 
p-{k) : 

p-{k)=f^{k) + -p1{k) + ... (4) 

In 3D PQ{k) is a volume term of order A;^, pf is a surface term of order k and the next order is 
a constant (order zero in k) edge term p2 • ^^ the following we specialize the periodic orbit theory 
initiated in Ref. ^ to the simple case of zero length trajectories and treat the boundary conditions by 
using the Balian and Bloch method [^. We will derive the formulae in the case of three dimensional 
billiards and only state the two dimensional results. 



The projector onto an irreducible invariant subspace (a) is written as [13| 



P'' = ^T.x''{9)uh9), (5) 

where |^| is the order of the group, x^id) is the character of group element g in irrep (a) and U{g) 
is the operator representing the action of g in Hilbert space. The projected level density is then (see 

Ik 

p^{k) = Im Tr{P"G(fc + iO+)} , (6) 

TT 

where G{k + iO"^) is the retarded Green function of our problem ; for simplicity we will drop the 
argument in the following. According to the ideas of Balian and Bloch, G is written as the free Green 
function Gq plus a correcting term : 

G = Go + Gi . (7) 

We will here only need the explicit expression of Gi for short trajectories near the boundary dB. 
In this case it can be approximated using the method of images : 

<r|Gi|r' >=e<r|Go|r'i > , (8) 

where e = ±1 for Dirichlet or Neumann boundary conditions and r'l is the symmetric of point r' with 
respect to the plane tangent to dB. Note that (^) is properly defined only if r and r' are close to 
each other and close to the boundary, (pi) is the leading correction (in the k expansion) to Gq arising 
near the boundary. A more systematic expansion can be found in Ref. |5|. In the Balian-Bloch 
terminology, we only consider the first term of the multiple scattering expansion and furthermore 
restrict ourselves to zero length orbits (or high wave vector). 

In three dimensions, Gq contributes to the total level density with a volume term /Oo(fc) = Vk"^ /Air'^ 
iy being the volume of the billiard), Gi gives the surface term pi{k) and the next order p2 is obtained 
by pursuing the multiple scattering expansion. We shall not consider higher orders in the expansion. 
Surprisingly enough, when computing the symmetry-projected terms (pg, pf and pf), we need not 
know the form of the Green function to a better approximation than (0,^). 

The contribution of Gq to the level density (|6|) is (using (^)) 

--TFiT. X'^ia) Im d\ < 9r\Go\r > , (9) 

TT \y\ '—' J 

I I g&g 

and the non-oscillating part of this expression is given by the quasi zero length trajectories going from 
r to gr. In 3D the elements of a point group Q can be either symmetries with respect to a plane, 
rotations or improper rotations (ie product of elements of the two first types) ; and orbits going from 
r to gr do not all contribute at the same order in the Weyl expansion : 

(i) If 5 = e is the identity, then er = r for all points in B and (^) gives the projected volume term 
-p%{k). 

{ii) If g is a symmetry with respect to a plane, it leaves invariant a surface of area Sg (intersection 
of the plane of symmetry with B) and then the corresponding term in (^ will contribute to /jf . 



(in) If 5( is a rotation, its invariant points in B are located on a segment of length Lg and we will 
have an edge contribution to P2 • 

Terms arising from improper rotations in (^ do not contribute at order pf s-iid "^^ '^ill ^ot consider 
them here. From (^ we can readily write the symmetry projected volume term (the first term in the 
Weyl expansion) : 

Poik) = ^-^Po(k) . (10) 

This formula is quite general and holds also for smooth potentials. A similar expression can be 
derived for continuous symmetry groups as well : the discrete sum in (|^ is simply replaced by an 
integration on the invariant measure of the group. To fix the ideas, if we consider a system with 5 
fold and 1 fold degenerate levels, formula (|l^) implies that roughly speaking, the 5 fold degenerate 
levels will contribute 25 times more to the level density than the singly degenerate levels (due to the 
coefficient (d")^ in (|To|)). This contribution comprises a trivial factor 5 coming from the degeneracy 
of the levels. But the factor 5 remaining (25 = 5 x 5) is more surprising : when choosing a level 
at random in the spectrum, the probability is approximatively 5 times bigger for drawing a 5 fold 
degenerate level than for a singly degenerate level. 

The contribution of Gi to the level density (|6[) is of the same type as Eq. (|9|) . There similarly, the 
identity, symmetries with respect to a plane and rotations do not contribute at the same order. The 
identity gives a surface contribution {d°')'^pi/\Q\ (as for the volume term (p!o|) arising from g = e in 
(^)). Symmetries contribute to P2 and we will not consider the effects of rotations (contributing at 
higher order). The term following Gi as given by (^) in the multiple reflexion expansion contributes 
to the total level density at order p2- Its contribution will be weighted by a factor (d")^/|^| when 
projected onto a given irrep (this again comes from gf = e in (||) while other elements of the group 
contribute at higher orders). 

Gathering all the contributions we get for the surface and edge terms : 

Pnk)=^-^Pi{k) ^ E x"(5)Im/dV<5r|Go|r>, (11) 

\y\ ^ t/ ^n J 



gGQsym. 



P2{k) = \^-P2{k) -^ Y. x"(5)Im/dV<5r|Go|r> 

\y\ TT g _^^^ J 



g€Grot. 

2k d'^ r 

- -Yg\ E x"(5)ImydV<5r|Gi|r> . (12) 

gdiQsym. 

In (11,12) the integration must be restricted to quasi-zero length trajectories. X^gGgrot. indicates a 
sum restricted to the rotations of the group (the same convention holds for X^gegsym.)- To be specific 
we will compute one of the integrals in (Il,|l2). Let us consider for instance the last term of the 
r.h.s. of Eq. (0). The integral /d^r can be separated in a surface integral Jdxdy along the invariant 
symmetry plane plus a term Jdz along the normal to this plane. Then the modulus \\gr — r|| = 2\z\ 
and the free Green function reads (see e.g. Ref. [|l4[ ) 



^ik\\gr-r\\ ^2ik\z\ 

<gr\Go\r>=--— =-^^. (13) 

47r||5rr-r|| 87r|2;| 

Hence the integral over dxdy will give a factor Sg (the area of the invariant surface) and the 
imaginary part of the integral over dz will give 

1 r-,/in(2fc|.|)^_l_ ^^^^ 



8tt J^oo k, 

Thus the total surface projected term ( |ll] ) reads : 

^""(^) = ^^"i(^) + S ^ X"(9)^. (15) 

gGysym. 

We note that : (i) if a class of G contains one symmetry it is only formed by symmetries ; (ii) 
X"ig)Sg is a class function, i.e. it takes the same value for all symmetries belonging to the same class 



C (we note x"(C)'S'c the value of the class function). Hence (15) can be written as a sum over those 



classes of the group which are formed by symmetries with respect to a plane (noted J2csym. below) : 

p?ik) = ^pi(fc) + ^ E x'-iQm^ . (16) 

\y\ \y\ ^ 47r 

In ( [l^ ) \Qc\ is the number of elements of the group belonging to class C. 
Similar computations yield for the edge terms : 

-P^k) = ^p-2(fc) + ^ E x"(C)|ac| ^ + ^^ E x"(C)|g.|f^ . (17) 

1^1 l^lcroi. 47rsin2(^) \^\csym. ^^ 

In (|l7| ) /3c is the angle of rotation, common (modulo an irrelevant sign) to all the rotations belonging 
to a given class C and Lc is the common length of the invariant segments. Here also we use the fact 
that if a class contains one rotation, then it contains only rotations. dSc is the common perimeter of 
the surfaces invariant under the symmetries belonging to the same class (i.e. dSc is the perimeter of 

Sc)- 

Formulae ([To|), (p^) and (|l^) (together with their 2D counterparts ([l8|) and (|l9|)) are the most 
important results of this paper. They allow to compute the first terms of the symmetry projected 
Weyl expansion knowing only the character table of the group and the surfaces, lengths and perimeters 
of sub-parts of B invariant under a given class of group transformations. Note that using elementary 



group theoretical properties of the character table (see e.g. Ref. |]13[) one recovers the total smooth 
level density when summing the different projected terms. 

In two dimensions, the present formalism allows to determine the value of Po{k) (which is now a 
surface term) and of Pi{k) (the perimeter correction). The formulae read : 

Poik) = ^Poik), (18) 

\y\ 

-pUk) = ^pi(^) + S E x'-mGci^. (19) 

\y\ \y\ r< 2^ 

' ' ' ' Csym. 



In (|T^|T^) po{k) = Sk/2'K {S being the area of the bihiard) is the total surface term and pi{k) = 
€L/4tt (L is the perimeter) the total perimeter term. Lc is the common length of the segments invariant 
under the symmetries belonging to the class considered. Note that similar results were obtained in 



Ref. 1 15] for a smooth potential (the two-dimensional quartic oscillator) using the Wigner transform 
approach. 

An alternative method for deriving the above results would be to work in reduce configuration 
space, to find out for each irrep which are the boundary conditions on the elementary cell and to work 



out the corresponding Weyl expansion (as was done in |1(:] for instance). But for some symmetry 
groups the boundary conditions on the elementary cell relevant for each irrep can be cumbersome 
(see below the example of the icosahedron) and also the corresponding Weyl expansion might not 
be known up to order pf (^) • some irreps lead to mixed Neumann-Dirichlet conditions or impose 
phase-shifts when going from one boundary of the elementary cell to an other (this is the case for the 
rather simple group C^v)- The method presented here has the advantage of giving simple and easily 
applied formulae that do not require detailed knowledge of the reduced boundary condition. 

As an illustration let us consider a three dimensional icosahedral billiard with Dirichlet boundary 
conditions. This billiard was studied in [|l^ as a model for faceted metal clusters. The first 565 
quantum levels were determined numerically (counting the degeneracies this leads to 2094 available 
states). The total symmetry group of the icosahedron is Ih ! to simplify the presentation we will 
consider here only its subgroup I formed by 60 rotations. X has 5 classes and 5 non equivalent irreps, 
we will present only the results for one of the irreps which has dimension d" = 3 (see table 1). In 
order to compare the numerical results with the predictions of Eq. ([To|), ( p!6| ) and ( p^ ) we focus on 
the cumulated level density : 

iV"(A;) = / p''{k')dk' = N^{k) + N^{k) + N^{k) + iV^(fc) + • • • (20) 

Jo 

From Ref. Q we know the total — i.e. summed over all the irreps — terms iVg, Ni and iV2 (the 
explicit expression is given in ||l^). I contains no symmetry and we can apply formulae ( [T^JT7| ) 
if knowing the length of the segments invariant under rotations (they are listed on table 1). The 
comparison of the numerical N°' — which is a staircase function — with its smooth approximation 



Nq + iVf + N2 (from (1C,16,13)) is shown Fig. 1. The difference is also plotted on this figure : as 
expected it is of order zero in k, and the amplitude of the oscillations is of order of the degeneracy 
d" = 3. A simple fit gives the value of the constant next order N^(k) ~ 0.288. The same agreement 
is also obtained for the other irreps of I or when studying the irreps of the total group X^. It is 
interesting to note that the boundary conditions to be applied on the elementary cell when studying 
the different irreps of X and Ih are very difficult to determine. The method exposed here uses the 
character table of the group which comprises the necessary information in an easily accessible form. 
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Table captions 



Table 1. Characters of the classes of T for the irrep considered in the text. For the designation 



of the classes we use the conventions of Ref. [18|. Lc is for each class (except for the identity) the 
length of the invariant segments. Lc is given in units of the edge length. 
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v^(3 - V5) 


y^Ws 


yiwi 



Figure captions 



Figure 1. bottom part : N°'{k) (staircase function) and its smooth approximation NQ+N1+N2 
as a function of k. k is given in units of the inverse length of the edge of the icosahedron. Top part 
: difference of the two quantities AN" = iV" - (iV^ + iVf + N^) as a function of iV°. The thin 
horizontal line is the mean value of AA^*^ (~ 0.288). The bottom plot is limited to the first 30 levels 
for legibility. The upper part concerns the first 100 levels {i.e. the first 300 quantum states). 



